Anomalous Hall effect in granular ferromagnetic metals 
and effects of weak localization 
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We theoretically investigate the anomalous Hall effect in a system of dense-packed ferromagnetic 
grains in the metallic regime. Using the formalism recently developed for the conventional Hall effect 
in granular metals, we calculate the residual anomalous Hall conductivity cr xy and resistivity p xy 
and weak localization corrections to them for both skew-scattering and side-jump mechanisms. We 
find that the scaling relation between p xy and the longitudinal resistivity p xx of the array does not 
hold, regardless of whether it is satisfied for the specific resistivities of the grain material or not. The 
weak localization corrections, however, are found to be in agreement with those for homogeneous 
metals. We discuss recent experimental data on the anomalous Hall effect in poly crystalline iron 
films in view of the obtained results. 

PACS numbers: 73.63.-b, 73.20.Fz, 61.46.Df 
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I. INTRODUCTION 

The anomalous Hall effect (AHE) in ferromagnetic ma- 
terials has been attracting the interest of researchers for 
decades. The first theoretical explanation 1 of AHE was 
given by Karplus and Luttinger in 1954. They have 
shown that, in essence, the anomalous Hall current arises 
from the population imbalance of the electron spin states 
that is transferred into the asymmetry in electron motion 
via spin-orbit coupling. Since then, the theory of AHE 
has undergone further significant developments (see, e.g., 
recent reviewsdand[^and references therein). The inter- 
pretation^ of AHE in terms of the Berry phase concept 
has fueled additional interesls^S^ to the problem. 

One distinguishes between the intrinsic and extrinsic 
AHE. The intrinsic AHE arises in a perfect periodic lat- 
tice subject to spin-orbit coupling, ft is due to the topo- 
logical properties of the Bloch states and does not require 
any disorder. On the contrary, the extrinsic AHE is due 
to the asymmetric spin-orbit scattering of spin-polarized 
electrons on the impurities of the sample. Two mech- 
anisms termed skew-scattering^ and side-jump^ are re- 
sponsible for the extrinsic AHE. They depend differently 
on the amount of disorder in the sample and, as a result, 
for certain type of disorder the anomalous Hall (AH) re- 
sistivity scales linearly (p xy oc p xx ) with the longitudi- 
nal resistivity for the skew-scattering, and quadratically 
(p xy ex p xx ) for the side-jump mechanism. These scaling 
relations were observed experimentally 1 ^ in homogeneous 
systems. At the same time, for some heterostructure sys- 
tem s 11 : 12 considerable deviations from this scaling law 
were reported. 

At sufficiently low temperatures, the physics of AHE 
is enriched by the quantum effects of Coulomb interac- 
tions and weak localization. The Coulomb interaction 
correction to the AH conductivity has been shown to 
vanish for both skew-scattering and side-jump mecha- 
nism s 13 : 14 . Weak localization (WL) effects were studied 
m Refs. [IHlIi and it was demonstrated that WL cor- 



rection to AH conductivity is nonzero 13 : 14 : 15 for skew- 
scattering and vanishe s] 14 : 15 for side-jump mechanism. 
The logarithmic temperature dependence of the AH resis- 
tivity and the absence of such for the AH conductivity ob- 
served in amorphous iron films 1 ^ were initially attributed 
to the Coulomb interactions^ 3 - and later interpreted^ 5 - in 
terms of the WL corrections for the side-jump mecha- 
nism. 

In a recent paper—, the logarithmic temperature de- 
pendence of the longitudinal and AH resistivities of the 
polycrystalline iron films at sufficiently low temperatures 
was reported. For well-conducting samples, the behav- 
ior could be well explained by the WL theor y 13 : 15 of the 
AHE in two-dimensional homogeneously disordered sam- 
ples. For more resistive samples, however, noticeable de- 
viations from the theoretical predictions were observed. 
The authors suggested that these deviations could be at- 
tributed to the granular structure of the samples. 

Motivated by the experimental data of Ref. [13, in the 
present paper we investigate AHE in a granular system 
of ferromagnetic metallic nanoparticles within a micro- 
scopic theory. For that purpose, we extended the recently 
developed theor y 18 : 19 of the conventional Hall effect in 
granular metals to describe AHE. 

The paper is organized as follows. In Sec. UH we for- 
mulate the model for a granular system. In Sec. IIII1 the 
residual AH resistivity is calculated, first using the clas- 
sical approach, and then this result is recovered from the 
diagrammatic approach. The breakdown of the scaling 
relation between the AH and longitudinal resistivities is 
discussed. In Sec. lIVl we calculate the WL corrections to 
the AH resistivity and discuss the experiment of Ref. IT?! 
Concluding remarks are presented in Sec. [V] 



II. MODEL 

We consider a regular quadratic (d = 2, single gran- 
ular monolayer) or cubic [d = 3, many monolayers) lat- 
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FIG. 1: Granular system and the classical picture of the Hall 
conductivity. The external Ohmic voltage V y is applied to the 
contacts in the y direction. The resulting Ohmic current I y = 
GrVy running through the grain in the y direction causes the 
Hall voltage drop Vb = RhI v between its opposite banks in 
the x direction. Since when calculating the Hall conductivity 
a xy the total voltage drop per lattice period in the x direction 
is assumed 0, the Hall voltage Vn is applied with an opposite 
sign to the contacts in the x direction (see bottom), causing 
the Hall current I x = G T V H = G\R u V y [see Eq. (pTTj) ]. 



tice of identical in form and size three-dimensional metal- 
lic grains coupled to each other by tunnel contacts with 
identical conductances Gt- At the same time, we as- 
sume that the grains are disordered either due to impu- 
rities in the bulk of the grains or due to an atomically 
irregular shape. The assumptions of the regularity of 
the system simplify the analysis significantly, but are not 
crucial. The results we obtain are expected to apply to 
structurally disordered granular arrays as well. We con- 
sider the metallic regime in this paper, when the tunnel 
conductance Gt = (2e 2 /H)gT is much larger than the 
quantum conductance, 



9T > 1. 



(2.1) 



In this limit, the whole granular system is a good con- 
ductor and the quantum effects of weak localization and 
Coulomb interactions can be studied 20 pcrturbativcly in 
1/<?t < 1. As usual, it is assumed that the granularity is 
well-pronounceo™, i.e., the dimensionless grain conduc- 
tance go exceeds the tunnel conductance gx, 



9o > 9t- 



(2.2) 



The key ingredients of AHE that give rise to a finite 
transversal conductivity a xy are (i) the spin magnetiza- 
tion of conduction electrons and (ii) considerable spin- 
orbit interaction. Analogously to homogenously disor- 
dered metals, the simplest Hamiltonian containing these 
two ingredients and thus describing AHE in a granular 
system can be written as 

H = H Q + Hu + H T . (2.3) 
The first two terms in Eq. (|2.3[) describe isolated grains, 



where 
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(2.4a) 

contains the kinetic energy and the exchange field h = 
he z directed along the z axis in all the grains (we put H = 
1). The exchange field causes a finite spin magnetization 
of the electrons. The intragrain disorder is described by 

Hu = Y.J dri¥t(rO{tf(n) 

i 

-iA 2 (a ■ [VUfa) x V n ])} *(n), (2.4b) 

where the first term corresponds to the conventional scat- 
tering on the disorder potential J7( r and the second one 
to the spin-orbit scattering. In Eqs. (|2.4ap and (|2.4b[) . 
v & = fai 7 ipi) T is the two-component spinor field oper- 
ator of the electrons, <r = (cr x , a y , a z ) denotes the vec- 
tor consisting of the Pauli matrices a a , a = x,y, z, and 
i = (ii, . . . , id) is an integer tuple numerating the grains 
on the lattice. The integration with respect to ri is per- 
formed over volume of the grain i. 

We consider the simplest model of disorder 



U(r) = ^2uS(r - r Q ), 



(2.4c) 



in which the point impurities are located at random posi- 
tions r a within the grains and are uniformly distributed 
with the concentration n\ over the volume of the grains. 
We assume that spin-orbit coupling is weak in the sense 
App <C 1, where is the Fermi momentum, and that the 
exchange field h ep is smaller than the Fermi energy 
€f of electrons in the grains. The latter two assump- 
tions allow one to study AHE perturbatively in h and 
spin-orbit coupling. 

The last term in the Hamiltonian (|2.3|) describes tun- 
neling between the grains, 

<iJ> 

Xy = Jdsid^tisus^fa)^). (2.4d) 

The summation in Eq. ()2.4d|) is done over the neighbor- 
ing grains i and j, so that each tunnel contact is counted 
only once, and the integration with respect to Si and sj 
is performed over two surfaces of the contact between 
the grains i and j, one belonging to grain i and the 
other to grain j. It is both physically reasonable and 
convenient for calculations^ to treat the tunneling am- 
plitudes i(si,Sj) as Gaussian random variables with the 
variance (tfa, Sj)t(sj, Si)) t = t§ 5(si — Sj). 

The anomalous Hall conductivity of the array is calcu- 
lated using the Kubo formula for granular systems in the 
Matsubara representation, 



CabM 



2-d 
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CO 



[IlabM - n ab(0)] , (2.5) 
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where 

nabM^-V/ 7 dre iwT (T T [7i,a(T)/j,b(0)]) (2.6) 
. Jo 

is the correlation function of the tunnel currents 

7 1 , a = -ie(X 1+a ,i-X M+a ). (2.7) 

Here, u> is a bosonic Matsubara frequency (we assume 
lu > throughout the paper), the lattice unit vectors a 
and b denote the directions of the current and external 
electric field, respectively. The approach to calculating 
the AH conductivity is analogous to that developed for 
the ordinary Hall effect in Ref. 0. It is based on the di- 
agrammatic perturbation theory in the tunnel Hamilto- 
nian (|2.4d[) with the ratio gr/go [Eq. (|2.2p ] of the tunnel 
and grain conductances as an expansion parameter. We 
refer the reader to Ref. for the details of the approach. 

In our model (l2~3l) - ([2~4dl) . the source of spin-orbit scat- 
tering are the impurities in the bulk of the grains [sec- 
ond term in Eq. (|2.4b[) ]. The anomalous Hall current, 
therefore, arises from the bulk of the grains. To perform 
explicit calculations, we will assume the intragrain dy- 
namics is diffusive, i.e., the mean free path I <C a is much 
smaller than the grain size a. In the opposite case of clean 
grains, surface scattering is dominant and spin-orbit scat- 
tering off the grain boundary could be the major source 
of AHE. If the boundary roughness can effectively be 
modeled by scattering on impurities in the bulk of the 
grains, our results may also be applicable to the arrays 
of ballistic grains with chaotic intragrain dynamics. 



III. RESIDUAL ANOMALOUS HALL 
RESISTIVITY 

A. Classical approach 

We start by calculating the residual anomalous Hall 
conductivity a xy and resistivity p xy of a granular ar- 
ray, neglecting quantum effects of weak localization and 
Coulomb interactions. Actually, as we show in this sub- 
section, as long as quantum effects are neglected the AH 
conductivity can be obtained by means of the classical 
electrodynamics without using the Kubo formula. The 
diagrammatic approach that will further allow us to in- 
clude quantum effects is presented in Sec. MI CI 

Within the classical approach, the granular array can 
be considered as a resistor network with the tunnel con- 
tacts viewed as surface resistors with conductance Gt- 
The AHE occurs inside the grains and is fully character- 
ized by the AH resistance Rh of the each grain. Given 
Rh, in the leading order in gr/go -C 1, one can easily 
arrive (Fig. []} at the expression 

Oxy = a 2 d G T Rn (3-1) 



for the residual AH conductivity of the granular array. 
Since the longitudinal conductivity equals 

& xx — a 2 ~ d Gr, (3.2) 
for the AH resistivity of the granular system we obtain 

p xy = ^L = a d - 2 R H . (3.3) 

^ XX 

The AH resistivity is, therefore, expressed solely through 
the Hall resistance Rh of a single grain and is indepen- 
dent of the tunnel conductance Gt, which determines the 
longitudinal resistivity [Eq. (|3.2p j. 

To get a further insight into the problem, one should 
specify Rh more explicitly. The electron transport in 
the diffusive grains can fully be described by the specific 
longitudinal af x and AH erf ^ conductivities of the grain 
material. The AH conductivity 

< = ^ SS + ( 3 - 4 ) 

is a sum of two contributions due to skew-scattering (ss) 
and side-jump (sj) mechanisms. Given <jf x ana ; a s^ one 
can find the anomalous Hall resistance Rh of the grain 
by solving the electrodynamics problem for the distri- 
bution of the electric potential in the grain 19 . Ana- 
lyzing this problem, one obtains that Rh is expressed 
through the specific AH resistivity of the grain material 
Pf y = v¥ y /(°xx) 2 in thc following way 

Rh = A H pf y /a, (3.5) 

where the numerical factor Ah < 1 is determined by 
the shape of the grain only. For simple grain geometries 
Ah can be found explicitly, e.g., Ah = 1 for cubic and 
Ah = tt/4 for spherical grains. 

As follows from Eqs. (|3.3[) and (|3.5p . the AH resistivity 
of a three-dimensional granular array (d — 3, many grain 
monolayers) 

Pxy = A H p% (3.6) 

is determined by the AH resistivity of the grain material, 
up to a geometrical numerical factor determined by the 
shape of the grain. The AH resistivity of a granular film 
(d = 2, one to several monolayers) 

Pxy = A H p%/d z (3.7) 

is obtained by dividing the 3D result (|3.6[) by the thick- 
ness d z of the film. 

The results (I3.1l)-(|3.7p are actually analogous to those 
obtained in Refs. Il8laiid [l9l for the conventional Hall 
effect in granular metals, with the AH resistivity of the 
grain material pfy entering the equations instead of the 
conventional Hall resistivity. Specifics of AHE is reflected 
in, e.g., the breakdown of the scaling relation between the 
AH and longitudinal resistivities, as discussed in the next 
subsection. 
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B. Breakdown of the scaling relation 

In homogeneously disordered systems, for certain types 
of disorder the AH and longitudinal resistivities obey the 
scaling relation 



Pxy oc p 



(3.8) 



with the exponent 7 = 1 for skew-scattering and 7 = 2 
for side-jump mechanisms. This scaling originates from 
the fact that spin-orbit scattering, which results in the 
transversal current, is caused by the same impurity po- 
tential U{r) [Eq. (|2~4b|) ]. scattering off which is respon- 
sible for the finite longitudinal resistivity. 

The scaling relation (|3 . 8|) holds for the model of identi- 
cal randomly placed short-range impurities, described by 
Eqs. (|2.4b|) and (f2T4c|) . Within this model, the longitu- 
dinal and AH conductivities of the grain material equal 



(3.9) 



37T7li1t 2 ' 



7 gr,ss _ 
' xy 



5tt 



CF 



(3.10) 



(3.11) 



Here, v is the density of states at the Fermi level for 
h = 0. As seen from Eqs. (|3.9|) - (|3.1ip , as the impurity 
concentration m is varied, the resistivities pf x = l/er|^. 
and pfy indeed change according to Eq. f|3 .8[) [it is im- 
plied in Eq. (|3.8p that the variation of p x r y and p x T x is 
caused by the change of n- u i.e., the amount of disorder, 
whereas the strength u of the scattering potential of sin- 
gle impurities is fixed]. The scaling relation (|3.8|) thus 
holds for the specific resistivities of the grain material. 
Although Eqs. (|3.9p - (|3.1ip are obtained for weak impu- 
rity scattering (Born approximation), it can be shown^ 
that the scaling law (|3.8[) holds for strong scattering as 
well, since the dependence on the impurity concentration 
remains the same. However, for more complicated type 
of disorder with stronger finite-range correlations of the 
disorder potential the scaling relation may be violated. 

Comparing Eqs. (|3.2p and (|3.3|) . we see that no scaling 
relation similar to (|3.8p between the AH p xy and longi- 
tudinal 



,d-2 



/G q 



(3.12) 



resistivities of the whole granular array holds. This result 
is actually not surprising, since the longitudinal and AH 
transport in granular systems are governed by different 
mechanisms: the former is due to tunneling through the 
potential barriers between the grains, whereas the latter 
is caused by spin-orbit scattering inside the grains. Ac- 
cording to Eqs. (|3.2p and (|3.3p . if the granularity is indeed 



pronounced [Eq. (|2.2|) ]. the AH resistivity should not vary 
much for samples with noticeably different longitudinal 
resistivities. One could say that for granular systems, 
the scaling relation (|3.8p with the exponent 7 = holds, 
independent of the dominant mechanism of AHE. In this 
context, we note that considerable deviations from the 
scaling law (|3.8p have previously been observed experi- 
mentally in several types of heterostructure system o 11 ' 12 , 
in which the longitudinal resistivity was also governed 
by the structural disorder (such as transparency of the 
interfaces) rather than by the intrinsic disorder of the 
ferromagnetic material. 



C. Residual anomalous Hall resistivity via 
diagrammatic approach 

The classical approach allows one to easily obtain 
Eq. (|3.ip for the residual AH conductivity and make some 
interesting conclusions about AH transport in granular 
metals at high enough temperatures. However, it has 
nothing to say about quantum effects of weak localization 
and Coulomb interactions, which set in at sufficiently low 
temperatures. To study these effects on the AH trans- 
port, a more sophisticated diagrammatic approach based 
on the Kubo formula (|2.5p is needed. Before we proceed 
with the quantum effects in Sec. lIVi we first demonstrate 
here how the classical result (|3.1I) is reproduced within 
the diagrammatic approach. 

As demonstrated in Ref. the key object of the di- 
agrammatic approach to the Hall effect in granular sys- 
tems is the intragrain diffuson, i.e., the two-particle elec- 
tron propagator of an isolated grain. It contains all the 
information about the specific mechanism of the Hall ef- 
fect. As usual, the diffuson is formally defined as the 
disorder-averaged product of two Greens's functions. In 
the presence of the exchange field and spin-orbit scat- 
tering the electron Green's functions are matrices in the 
spin space and the intragrain diffuson is defined as their 
direct product, 

D u (r,r') = -L(4 +W (r,r')®^(r',r))r/, ( £ + u ) e < 

(3.13) 

Here, C/'s are the exact Green's functions of the intragrain 
Hamiltonian Hq + Hjj in the Matsubara technique for a 
given realization of the disorder potential U (r) and the 
angle brackets (...)[/ denote disorder-averaging. 

According to the Kubo formula (|2.5p . the conductivity 
is in the leading order expressed through the spin-singlet 
diffuson component, which is given by the trace Tr CT of 
the Green's functions in the spin space 

D u (r,r') = J-(Tv a [G £+ Ur,r')G £ (r' ,r)])u, {e+uj)e < 0. 

(3.14) 

Below, we will need the spin-singlet diffuson (|3.14p only. 

Analogously to the conventional Hall effec t 18 ' 19 , the 
residual AH conductivity a xy is given by the diagrams 
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FIG. 2: Diagrams for the residual anomalous Hall conductiv- 
ity a xy [Eqs. (EH! and (335)], see Ref. QJ for details. 

in Fig. [2] Calculating these diagrams, one can relate the 
tunnel conductance in Eq. (|3.1[) for <j xy to the microscopic 
parameters of the model as gr — 27r(z^o) 2 So (So is the 
area of the contact) and express the AH resistance of the 
grain through to the intragrain diffuson (|3.14p as 

Rh = ^(Ds-D^ + D/-D^), (3.15) 

where 

D <* = "^2 /dsods a £) w=0 (so,s a ) (3.16) 

are the diffusons at zero frequency u> = connecting 
different contacts as shown in Fig. [21 with a = 1,2,3,4 
for a =/*, \, y, \, respectively. 

The problem of calculating i?# is, therefore, reduced 
to finding the diffuson. Within the conventional disorder- 
averaging technique^, the diffuson (|3 . 14[) can be shown 
to satisfy the diffusion equation 

(w-A,V;)D u (r,r / ) = a(r-r / ), ( 3 - 17 ) 

in which Dq = v f t/3 is the coefficient of the intragrain 
diffusion (vp is the Fermi velocity and r is the scattering 
time, 1/r = 2nvriiU 2 ). Equation (|3.17p itself clearly does 
not contain any information about the Hall effect. This 
information is contained in the boundary condition for 
D u (r, r'), which Eq. (13 . 1 T|) must be supplied with for a 
finite system. In Ref. a general method of deriving the 
boundary condition for the diffuson was developed and it 
was shown that the boundary condition may be written 
as 

n a (j a rp)V r0 D u (r,r%es = 0. (3.18) 

Here, a,j3 = x,y,z, and n a are the components of the 
unit vector n normal to the grain boundary S at point r 




FIG. 3: Diagrams for the Hall component {jxr y ) BB of the 
current-coordinate correlation function p. 191) due to the 
skew-scattering mechanism. 

and pointing out of the grain. In Eq. (|3.18[) . 

(j a rp) = ^Tr CT J dr' (j ra [Qe+U^ r')S e (r', r)](r' - r)^ 

is the current-coordinate correlation function. The non- 
relativistic part of the current operator j has the conven- 
tional form 

j[^*(r),V(r)] = -^-[V*(r)W(r) - ^(r)V^*(r)] (3.20) 
Zm 

Explicit form of the boundary condition (|3.18[) is thus 
determined solely by (j a r^). Analogously to the con- 
ductivity tensor, only the longitudinal (j x T x ) — (jyfy) — 
(j z r z ) and Hall {j x fy) — ~(jy r x) components are nonzero 
(we remind the reader that the exchange field is directed 
along the z axis) . This allows us to rewrite Eq. (|3.18p in 
the form 

(j*ra)(n-V r ZUr,r'))| r€ s = (j x r y )(t ■ V r D w (r, r'))|res, 

(3.21) 

where the vector t = [n x h] /h is tangent the grain bound- 
ary at point r. 

Since the AHE is weak due to the smallness of the 
spin-orbit coupling constant \pp -C 1 and the exchange 
field hjtp <C 1, the longitudinal component (j x r x ) can 
be calculated neglecting the exchange field and spin-orbit 
scattering completely and the expression for it reads 

(j x r x ) = ~y (3 ' 22) 

All specifics of the AHE is contained in the Hall com- 
ponent (j x r y ). Analogously to the AH conductivity of 
a homogeneously disordered metal [see Eqs. (|3.4p . (|3.10p 
and (|3.1ip ] , the total Hall correlation function 

(j x r y ) = ( Jx r y ) ss + {j x r v yi (3.23) 

is the sum of two contributions due to skew-scattering 
(ss) and side-jump (sj) mechanisms. 

The skew-scattering part (j x f y ) is given by the di- 
agrams in Fig. [3l which contain the impurity lines de- 
scribing the third-order scattering processes on a single 
impurity, see, e.g., Ref. [1. Calculating these diagrams, 
we obtain 

(jxr y ) ss = ~^-vu(Xp F ) 2 — (j x r x ) . (3.24) 
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(a) 



(b) 



l r /3 





Having established the correspondence between the 
classical and diagrammatic approaches, comparing 
Eqs. ([535]) . (j3~T7j) . and (f^Tj) with Eq. ([311) and using 
the Einstein relation a& T x — 2e 2 vDo, we can express the 
Hall resistance of the grain as 



R H = A 



H 



1 UxTy) 1 

2e 2 v (j x r x ) D a 



(3.30) 



This form will be used in the next section for calculating 
WL corrections. 




FIG. 4: Diagrams for Hall component {jxr y ) sl of the current- 
coordinate correlation function (|3.19|) due to the side-jump 
mechanism. Diagrams (a) and (b) contain the relativistic 
corrections to the current [Sj, Eq. (|3,25[) ] and coordinate [<5r , 
Eq. (I3.26P ] vertices, respectively. 



The diagrams for the side-jump contribution (j x r y ) sl 
are shown in Fig. 2] The diagram in Fig. 2] (a) contains 
the conventional for the side-jump mechanism relativistic 
correction 

6} = A 2 [a x V£/(r)] , (3.25) 

to the current operator (|3.20[) . see, e.g., Ref. 0. Addi- 
tionally, there exists an analogous relativistic correction 
to the coordinate vertex. This contribution can be ob- 
tained by repeating the derivation of the boundary con- 
dition (|3.18p done in Ref. |l£j, but taking into account the 
spin-orbit term of Hjj [Eq. (|2.4b|) ] in the diffuson ladder. 
This gives the diagram in Fig. [?] (b) , in which 



Si = A 2 [<x x p] 



(3.26) 



is the relativistic correction to the coordinate operator 
(p = — iV). One can recognize that Sr is the operator 
of the lateral translation ("side-jump"), see, e.g. Ref. I22I. 
Calculating the diagrams in Fig. [U we obtain 



(j x r y f = -3tt(H 2 — {\p F f- (j x 

U6F £F 



(3.27) 



Compari ng E qs. (pT22"lh pT24|) and (|3~27|) with 
Eqs. ([5T9")) , (|3~TU|) and (|3~TT|) . we note that for both skew- 
scattering and side-jump contributions the relation 



\ss/sj 



s/sj 



(j x r x ) 



(3.28) 



holds. Therefore the boundary condition (|3.18j) may be 
rewritten as 

af x {n ■ V T D„(r, r'))| reS - 0* (t • V r D u (r, v'))\ reS . 

(3.29) 

As shown in Ref. [ijj, it is exactly this form of the bound- 
ary condition, which is necessary to reproduce the clas- 
sical result (13. 5p for the Hall resistance Rh obtained by 
solving the electrodynamics problem. 



IV. WEAK LOCALIZATION CORRECTIONS 

A. Calculations 

We now proceed with calculating the weak localization 
corrections to the obtained "classical" anomalous Hall 
conductivity (|3.1[) and resistivity (|3.3|) of the granular 
metal. 

Technically, one has to consider WL corrections to the 
diagrams in Fig. [2] for the bare Hall conductivity (|3.1j) 
by inserting the Cooperon ladders into them in all possi- 
ble ways. As shown in Ref. UH, such WL corrections are 
factorized according to the form Eq. (|3.ip . i.e., there are 
diagrams describing the corrections to the tunnel con- 
ductance Gt only and to the Hall resistance Rh of the 
grain only. This allows one to write down the total weak 
localization correction Sa xy to AH conductivity a xy in 
the form 



8a x 



5Gt SRh 
> _i_ 

Gt Rh 



(4.1) 



Naturally, the WL correction 5Gt to the tunneling con- 
ductance has the same form as that to the longitudinal 
conductivity a xx [Eq. (|3.2p ] 23 ' 24 and reads 



SGj 
Gt 
Here. 

CM) 



6a x 



^■[C &;= o(i + a J i) + C w= o(i,i + a)] (4.2) 



a d d d q 



a ia(q-(i-j)) 



(2tt) 



+ 2r^Jl-cos( 9a a)] + l/T, 

(4.3) 

is the Cooperon of the whole granular array calculated in 
the zero-mode approximation for the intragrain Cooper- 
on o 23 ' 24 i 25 , A is the mean level spacing in the grain, and 
a = e T or a = ej,. In Eq. (|4.3|) . T — $tA is the tunneling 
rate, the dephasing time t v was introduced by hand, and 
the integration with respect to the quasimomentum q is 
performed over the first Brillouin zone q £ [—it /a, Tr/a] d 
of the grain lattice. In order not to complicate the anal- 
ysis, we assumed in Eq. (|4.3p that the dephasing rate 
1/ 7 V ^ l/ T so exceeds the spin orbit scattering rate l/r so . 
If 1 /r„, 1 /t so , and h are of the same order, the spin struc- 
ture of the Cooperon can be taken into account as, e.g., 
in Ref. [H. 
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FIG. 5: Diagrams for the weak localization correction S(j x r y ) BB to the skew-scattering Hall component (j x r y ) BB [Eq. (I3.24|l ] of 
the current-coordinate correlation function ()3.19p . The gray regions denote the Cooperons [Eq. (|4.3j) ]. Each diagram depicted 
above stands for four diagrams: the ones not shown are obtained by reflecting the impurity lines through the diagram center 
and/or flipping the diagram upside-down. 




FIG. 6: Diagrams for the weak localization correction 8{j x r y ) BB to the side-jump Hall component {jxr y ) B1 [Eq. (|3.27|l ] of the 
current-coordinate correlation function t|3. 19|1 . Diagrams (a) and (b) contain the relativistic corrections to the current [<5j, 
Eq. (I3.25P ] and coordinate [St, Eq. (|3,26[l ] vertices, respectively. Each diagram in panel (a) stands for two diagrams: the one 
not shown is obtained by flipping the diagram upside-down. 



According to Eq. (|3.30|) , the AH resistance of the 
grain has been expressed through the diffusion coef- 
ficient Dq and the correlation functions (jx^x) and 
(jxfy), which fully characterize the intragrain diffuson 
A^r, r') [Eqs. l|5TT7|> and IpTIS]) ]. As these three are 
well-defined correlation functions, one can calculate the 
WL corrections 5D , 5{j x r x ), and 5{j x r y ) to them using 
the diagrammatic technique. This will allow us to ob- 
tain WL correction Sp xy to the AH resistivity p xy from 
Eqs. Q and ([330)) as follows 



Spa 



8R 



H 



5{j x r y ) S(j x r x ) 5 D 



Pxy 



Rh 



{jxTy 



(jxr x ) 



Dn 



(4.4) 



The WL corrections to the diffusion constant Do 
and longitudinal current-coordinate correlation func- 



tion (j x r x ) are identical to those in Ref. [19| and have 
the form 



SDp _ S(j x r x ) 

D {jxTx) 



(4.5) 



Here. 



c = — C w =a(i,i) 

IT 

A f a d d d q 



,(4.6) 



7T J (2n) d 2TJ2 a [l - cos(q a a)] + 1/t, 

All specifics of AHE is contained in the Hall compo- 
nent (j x r y ). The diagrams describing the WL corrections 
to (j x r y ) ss [Eq. (J3JMJ)] and (j x r y ) si [Eq. (j3~771) ] are ob- 
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tained from the diagrams in Figs. [3] and [4] by inserting 
the Cooperon ladder into them in all possible ways. 



to the 

skew-scattering correlation function {jxfy) ss - The dia- 



Let us first consider WL correction 6{j 



x< yl 



grams for 5(j x 



are shown in Fig. [5] In total, there 



are 32 diagrams. After a tedious but straightforward cal- 
culation, we obtain 



(jxTyY 



-C. 



(4.7) 



The diagrams for the WL correction 5{j x r y ) s i to the 
side-jump correlation function {j x fy) Si are shown in 
Fig. [H The total number of these diagrams is 13. Calcu- 
lating them, we find that the contributions from all these 
diagrams cancel each other identically, which results in a 
vanishing correction 



S(jx 



0. 



(4.8) 



As seen from Eqs. (|4.7p and (|4.8p . the results for WL cor- 
rection differ for skew-scattering and side-jump mecha- 
nisms. Inserting Eqs. (|3~5]) . P~7|) . and P~g|) into Eq. P~4"]l . 
for the WL correction to the AH resistivity we obtain 



„ss 
r.r.y 



$j4y 
Pxy 



= 2c 



(4.9) 



(4.10) 



for skew-scattering and side-jump mechanisms, respec- 



tively. The total AH resistivity p xy 



P% is the 



sum of the skew-scattering and side-jump contributions 
and for the total WL correction Sp xy = 8p s xy + Sp^ y one 
obtains from Eqs. fO]) and (|4.10j) that 



Spxy 
Pxy 



AxyCj 



where 



A - 
-**-xy — 



/> ss -I- Iff 'i 
rxy 1 rxy 



rxy 



Pxy 



(4.11) 



(4.12) 



The factor (|4.12|) equals A xy — 1 and A xy = 2 for 
prevailing skew-scattering (p^ y ^> p s J y ) and side-jump 
(p x S y ^ Pxy) mechanisms, respectively, and belongs to 
the range 1 < A xy < 2, when two mechanisms give com- 
parable contributions. 

Equations (|4"TTj) and P~T2"]) . together with Eq. (jH)]) . 
constitute our final result for the weak localization cor- 
rections to the anomalous Hall resistivity of a granular 
metal. In the next subsection, we discuss the obtained 
results. 



B. Discussion 



We note that the form of Eqs. (|4~6]) , ([iJ]) . and (|4~10l) 
agrees with the results for WL corrections to the AH 
resistivity of homogeneous metal a 13 i 15 for both skew- 
scattering and side-jump mechanisms. This is most 
clearly seen, when the main contribution to the integral 
over q in Eq. (|4.6| comes from small momenta, qa <C 1, or 
equivalently, from spatial scales much exceeding the grain 
size. This happens in two (d = 2) and one {d = 1) dimen- 
sions, the latter case of granular "wires" is, however, irrel- 
evant for the Hall effect. In three dimensions (d — 3), the 
integral over q in Eq. (|4.6[) converges, if one neglects de- 
phasing, and the relative correction c(t v ) ~ 1/<7t, there- 
fore, depends only weakly on t v for finite dephasing. 

For two-dimensional arrays (one to several grain mono- 
layers), neglecting dephasing, the integral with respect to 
q in Eq. f|4. 6|) is logarithmically divergent at small mo- 
menta ga < 1. This divergency is cut by the finite de- 
phasing rate 1/t^. At low enough temperatures, when 
1/t v <C r, the divergency is strong and one obtains 



27r.9n 



ln(r-r v ), 



(4.13) 



where the dimensionless sheet conductance gu = 
axx/l^e 2 / (2irh)] — 2irgT of the array was introduced. 
At higher temperatures, when the dephasing rate \/t v 
becomes of order or larger than the tunneling rate T, 
1/t v > T, the integral in Eq. (|4.6|) is not divergent at 
qa <C 1 and the WL corrections of the granular film are 
not logarithmic in r„ anymore. At even higher temper- 
atures, when the dephasing rate exceeds the Thouless 
energy of the grain, 1/r^ > Do /a 2 , the contributions to 
WL corrections come from the bulk of each single grain, 
whereas the coherence of the intergrain motion is de- 
stroyed. Since in realistic granular systems the grains 
are three-dimensional particles, the WL corrections in 
this case are given by the result^ for a three-dimensional 
sample, 



5p X y 

Pxy 



— A xv 



_3\/3_ 

(e F r) 2 



1/t v > D /a 2 



(4.14) 

Here, a ~ 1 is a numerical cutoff-dependent factor and 
r is the intragrain scattering time. The correction (|4.14|) 
has a conventional for the 3D case square-root depen- 
dence on the dephasing rate. So, the "large-scale" low- 
temperature regime 1/t v <C T is the only one, in which 
the WL corrections of a granular film are logarithmic in 
t v . 

Using Eq. (|4. 13[) . we can write down the WL correction 
(l4"TTj) in the limit 1 /t v <C T as 



5p xy 

Pxy 



A, 



1h(Tt ¥ ). 



(4.15) 



In this form, the result (14. 1 5[) is in full agreement (up 
to a different infrared cutoff scale T, which is deter- 
mined by the microscopic structure of the system) with 
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that for a conventional homogeneously disordered metal 
characterized by the same sheet conductance g^. This 
sort of "universality" is actually quite expected, since 
WL corrections in 2D arise from large mesoscopic spatial 
scales, at which the microscopic structure of the mate- 
rial, whether it is homogeneous or granular, becomes ir- 
relevant. Therefore, it would be impossible to distinguish 
between granular and homogeneous two-dimensional ma- 
terial by measuring WL corrections. In this context, we 
remind that WL correction to the longitudinal 2 — and con- 
ventional Hall 19 resistivities of a granular metal have ear- 
lier been shown to agree with those for homogeneously 
disordered metals. In line with Eq. (|4. 15[) . one can write 
down WL correction to the longitudinal resistivity of a 
granular film in the form 2 ^ 



Pxx 



A x 



in(r v ). 



(4.16) 



with A xx = L 

In view of the obtained results, we would like to dis- 
cuss the recent experiment of Ref. [It]- The authors of 
Ref. [13 reported on the logarithmic temperature depen- 
dence of the longitudinal p xx (T) = R xx ln(Xb /T) and AH 
Pxy{T) = R xy \n(To/T) resistivities of the polycrystalline 
iron films at sufficiently low temperatures. For the most 
conductive samples, the values of the prefactors R xx and 
R xy were in a good agreement with the theoretical pre- 
diction o 13 ' 15 for WL corrections in two-dimensional ho- 
mogeneously disordered metals for the case of the dom- 
inant skew-scattering (provided one assumes the linear 
1/t v (T) cx T temperature dependence of the dephasing 
rate, as predicted, e.g., for electron-electron interactions 
by the diffusive Fermi liquid theory for both homoge- 
neous 2 ^ and granular— metals). This means that the 
In T-dependencies of p X x{T) and p xy {T) were well de- 
scribed by Eqs. (|4~T5f and (|4~16f with A xy = 2 (indi- 
cating that side-jump mechanism of AHE is dominant in 
these samples) and A xx = 1- This suggested the expla- 
nation of the observed behavior in terms of WL effects. 
For more resistive samples the In T-dependence seemed 
to persist, but the prefactors deviated significantly from 
the predicte d 13 i 15 values. That is, the behavior of p xx {T) 
and p xy (T) could still be described by Eqs. (|4. 15|) and 
(|4.16[) with, however, smaller prefactors A xy < 2 and 
A xx < 1 . The authors argued that these deviations could 
be explained by the onset of granularity in more resistive 
samples. 

According to Ref. [Tt|, in the regime of the intragrain 
dephasing length l v <C a {1/t v = Dq/1^) smaller than 
the grain size a, the WL correction to the AH resis- 
tivity p xy had to be given by Eq. (|4.15|) . but with 
the grain conductance go entering the denominator of 
the prefactor instead of the tunnel conductance gx, 
Sp X y/p xy oc A xy hx{ltp/l) I go. This would indeed be the 
case for flat pancake-shaped grains provided their 2D 
size a were much greater than their thickness ao, so that 
ao *C l v <C a. Considering that the most resistive sam- 
ples in the experiment of Ref. [Tt] were about 2nm thick, 



this would require the grain size a to be at least 20nm. 
However, the authors of Ref. [l?] presented an estimate 
of the tunneling rate T for lnm grains, which would cor- 
respond to the case of 3D grains. In such a situation, 
the WL correction to the AH resistivity should be de- 
scribed by Eq. (|4.14[) and it is not logarithmic in t v . 
Moreover, using the method developed in the present pa- 
per one can demonstrate that for pancake grains in the 
regime ao <C l v <C a the WL correction to the longitu- 
dinal conductivity would be a sum of the logarithm and 
logarithm squared contributions in the dephasing length, 
Sa xx /a xx oc -g T /3 2 [ln 2 (^/Z ) + (3\ri(l v /l)) (/3 ~ 1 is a 
geometrical factor, Iq is the contact size). This yields 
the In 2 T -dependence of the longitudinal resistivity p xx , 
which does not seem to agree with the data of Ref. [l7l . 
where both p xx and p xy were logarithmic in temperature. 
For these reasons, we do not think that the model of pan- 
cake grains corresponds to the experimental situation of 
Ref. M- 

At the same time, in the limit 1/t v 3> T, we have 
demonstrated for the AH resistivity [Eq. (|4. 1 5[1 ] and it 
has been earlier shown^ 3 - for the longitudinal resistiv- 
ity [Eq. (|4.16p ] that WL corrections are essentially the 
same for granular and homogeneously disordered metals. 
Since this is the only regime, in which the WL correc- 
tions to both pxx and p xy of a granular film are logarith- 
mic in t v , we conclude that the observed deviations of 
the prefactors A xy and A xx from the values A xy = 2 and 
A xx = 1 cannot be explained by the granular structure 
of the system and one should find an alternative expla- 
nation of the effect. We emphasize that in Ref. [I?] not 
only the prefactor A xy for the Hall, but also for the lon- 
gitudinal resistivity deviated from its "universal" value 
A xx = 1 . Since AHE in the experiment of Ref. [l?] is weak 
in the sense p xy <C p xx , the longitudinal resistivity p xx is 
not noticeably affected by the Hall effect. Therefore, the 
conclusion that the WL effects in a granular metal can- 
not explain the observed behavior could be drawn based 
alone on the earlier result (|4.16[) for the longitudinal re- 
sistivity, without any knowledge about AHE. 

Let us also briefly discuss the role of the Coulomb in- 
teractions in context of the data of Ref. [I?]- In Refs. 
and HO, the Coulomb interaction corrections to the lon- 
gitudinal resistivity that are specific to granular metals 
and absent in conventional metals were found. Analo- 
gous corrections were shown to exist for the conventional 
Hall resistivity in Ref. [l^ and one could demonstrate that 
the results of Ref. [li^ also apply to the AH resistivity. As 
these Coulomb interaction corrections are logarithmic in 
temperature (in any dimensionality of the array), one 
could be tempted to explain the data of Ref. [T3 in terms 
of them. Unfortunately, this would not be possible, since 
these corrections are of insulating nature, i.e., the relative 
corrections to the resistivities p xx and p xy are positive. 
Therefore, taking them into account would increase the 
value of the prefactors in the logarithmic T-dependencies 
of the AH and longitudinal resistivities. This would be in 
contradiction with the data of Ref. [13, where a decrease 
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of the prefactors A xy and A xx for more resistive samples 
was observed. 



V. CONCLUSION 

In conclusion, we have theoretically investigated the 
anomalous Hall effect in ferromagnetic granular metals. 
We found that no scaling law relation between the resid- 
ual anomalous Hall and longitudinal resistivities of a 
granular metal holds, regardless of whether this scal- 
ing holds for the specific resistivities of the grain ma- 



terial or not: the Hall resistivity of the whole array does 
not change as the longitudinal resistivity of the array is 
varied. At the same time, the weak localization correc- 
tions to the anomalous Hall resistivity of two-dimensional 
granular metals are found to be in full agreement with 
those for conventional metals. This is explained by the 
fact that the weak localization effects in low-dimensional 
conductors are determined by large mesoscopic spatial 
scales, at which the microscopic structure of the system 
is indistinguishable. 

Financial support of SFB Transregio 12 is greatly ap- 
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